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ABSTRACT 



Context. X-ray observations performed with the Rontgen Satellite (ROSAT) led to the discovery of seven radio-silent 
isolated neutron stars (INSs) which are detected only through the relatively dim and purely thermal X-ray emission 
from the cooling star surface. A few of these INSs (a.k.a. X-ray Dim INSs, or XDINSs) have been also detected at 
optical wavelengths where they seem to feature thermal spectra. Optical studies of XDINSs thus play a crucial role in 
mapping the temperature distribution on the neutron star surface and in investigating the existence of an atmosphere 
around the neutron star. 

Aims. The aim of this work is to investigate the optical identification of the XDINS RX J0420. 0—5022, tentatively 
proposed by Haberl et al. (2004) based on Very Large Telescope ( VLT) observations. 

Methods. We re-analysed the original observations of Haberl et al. (2004) to assess the detection significance of the 
proposed counterpart and we performed deeper VLT observations aiming at a higher confidence detection. 
Results. With a ~ 2a detection significance and a re-computed flux of B — 27.52 ± 0.61, we can not rule out that the 
proposed counterpart was spurious and produced by the halo of a very bright nearby star. While we could not detect the 
proposed counterpart in our deeper VLT observations, we found evidence for a marginally significant (~ 3.9cr) detection 
of a similarly faint object (B — 27.5 ± 0.3), w 0'.'5 north of it and coincident with the updated Chandra position of 
RX J0420. 0—5022. Interestingly, the angular separation is consistent with the upper limit on the RX J0420. 0—5022 
proper motion (Motch et al. 2009), which suggests that we might have actually detected the Haberl et al. proposed 
counterpart. From the flux of the putative RX J0420. 0—5022 counterpart we can rule out a > 7 optical excess with 
respect to the extrapolation of the XMM-Newton spectrum. 

Conclusions. High spatial resolution observations with the refurbished Hubble Space Telescope (HST) are the only way 
to confirm the detection of the putative candidate counterpart and to validate its identification with RX J0420. 0—5022. 

Key words. Optical: stars; neutron stars: individual RX J0420. 0—5022 



1. Introduction 

X-ray observations performed with the Rontgen Satellite 
(ROSAT) yielded to the identification of a group of seven 
radio-silent (Kondriatev et al. 2008j3 Isolated Neutron 
Stars (INSs). Their relatively dim X-ray emission (Lx ~ 
10 30 -10 31 ergs -1 ) originally earned them the nickname 
of X-ray Dim INSs, or XDINSs (see Haberl 2007; van 
Kerkwijk & Kaplan 2007, and Kaplan 2008 for recent re- 
views). Recently, a new XDINS candidate has been iden- 
tified in archival XMM-Newton observations (Pires et al. 
2009). XDINSs have purely thermal X-ray spectra which 
are best represented by a blackbody (kT w 50-100 eV), as 
expected for middle-aged (~ 1 Myr) cooling INSs, whose 
emission radius is consistent with a sizable fraction of the 
neutron star surface. The derived hydrogen column densi- 
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ties Nn ss 10 cm suggest distances < 500 pc (Posselt et 
al. 2007), as confirmed in two cases by their optical paral- 
laxes (e.g. van Kerkwijk & Kaplan 2007). X-ray pulsations 
(P = 3-12 s) have been detected for all of them (Haberl et 
al. 1997, 1999; Haberl & Zavlin 2002; Hambaryan et al. 
2002; Zane et al. 2005; Tiengo & Mereghetti 2007) but 
RX J1605. 3+3249, although with different pulsed fractions. 
The measurement of the period derivative P (Cropper et 
al. 2004; Kaplan & van Kerkwijk 2005a; Kaplan & van 
Kerkwijk 2005b; van Kerkwijk & Kaplan 2008) yielded 
spin-down ages of ~ 1.5 — 3.8 Myrs and rotational energy 
losses E ~ (3 — 5) x 10 30 erg s _1 . Broad absorption features 
(E\ ine 0.2-0.7 keV) have been observed in all XDINSs but 
RX J1856. 5-3754 (Haberl et al. 2003, 2004; van Kerkvijk 
et al. 2004; Zane et al. 2005), superimposed to the ther- 
mal continuum. These features are likely due to proton cy- 
clotron and/or bound-free, bound-bound transitions in H, 
H-like and He-like atoms. The inferred magnetic fields of 
~ 10 13 — 10 14 G are consistent with the values derived from 
the neutron star spin down and suggest that XDINSs might 
be (evolutionary) linked to other class of INSs, the magne- 
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tar candidates (see Mereghetti 2008 for a recent review) 
and the Rotating Radio Transients (e.g. Popov, Turolla, 
Possenti 2006). 

In the optical, only RX J1856.5-3754 (Walter & Matthews 
1997; Walter 2001), RX J0720.4-3125 (Motch & Haberl 
1998; Kulkarni & van Kerkwijk; 1998; Motch et al. 2003) 
and RX J1605. 3+3249 (Kaplan et al. 2003; Motch et al. 
2005; Zane et al. 2006) have counterparts certified by their 
proper motion measurements, while likely candidates have 
been proposed for RX J1308.6+2127 (Kaplan et al. 2002) 
and 1RXS J214303.7+065419 (Zane et al. 2008; Schwope 
et al. 2009) based on their coincidence with the X-ray posi- 
tions. Apart from providing a clear evidence of the optical 
identification, proper motion measurements are important 
to obtain an estimate of the kinematic age of the neutron 
star, to be compared with the characteristic age derived 
from the spin-down. The XDINS optical fluxes usually ex- 
ceed by a factor of ~ 5 (or more) the extrapolation of the 
X-ray blackbody, and their optical spectra, when measured, 
seem to follow a Ray leigh- Jeans distribution (e.g., Kaplan 
2008). The XDINS optical emission has been interpreted 
either in terms of a non-homegcneous surface tempera- 
ture distribution, with the cooler part emitting the optical 
(e.g., Pons et al. 2002), or of reprocessing of the surface 
radiation by a thin H atmosphere around a bare neutron 
star (Zane et al. 2004; Ho 2007), or of non-thermal emis- 
sion from particles in the star magnetosphere (Motch et al. 
2003). However, for the measured E, magnctospheric emis- 
sion would not be detectable, at least if an average optical 
emission efficiency of rotation-powered neutron stars (e.g. 
Zharikov et al. 2006) is assumed. Alternatively, like for the 
magnctars, optical magnetospheric emission might be pow- 
ered by the neutron star magnetic field, as proposed for 
1RXS J214303. 7+065419 (Zane et al. 2008). 
One of the XDINSs without a certified optical counterpart 
is RX J0420. 0-5022. The first optical observations of the 
field performed with the New Technology Telescope (NTT) 
soon after the discovery of the X-ray source (Haberl et al. 
1999) did not reveal any candidate counterpart brighter 
than 25.2 and R^ 25.2. More recently, thanks to the 
updated Chandra position, a possible optical identification 
was proposed by Haberl et al. (2004) with a faint object 
(B= 26.6 ± 0.3, V> 25.5) tentatively detected on archival 
Very Large Telescope ( VLT) images. However, the identi- 
fication has not been confirmed so far. The field of RX 
J0420. 0—5022 was also observed in the near-infrared (NIR) 
with the VLT but no candidate counterpart was detected 
down to H~ 21.7 (Mignani et al. 2007; Lo Curto et al. 2007; 
Posselt et al. 2009) and K s - 21.5 (Mignani et al. 2008). 
In this paper we re-analyze the original VLT observations of 
RX J0420.0-5022 presented by Haberl et al. (2004) and we 
report on follow-up, longer optical observations of the can- 
didate counterpart, performed by our team with the VLT. 
Observations and data analysis are described in Sect. 2, 
while results are presented and discussed in Sect. 3 and 
Sect. 4, respectively. 

2. Observations 

2.1. Observation description 

Optical observations of RX J0420. 0—5022 were performed 
in service mode with the VLT at the ESO Paranal obser- 
vatory on November 21st 2000, on November 25th 2006, 



Table 1. Log of the VLT FORS1 and FORS2 B-band ob- 
servations of RX J0420. 0—5022). Columns report the ob- 
serving date (yyyy-mm-dd), the number of exposures (N) 
and the total integration time per night (T), the image 
quality (IQ) and rms (in parentheses), as computed on the 
image, and the airmass. Values are the average computed 
over the exposure sequence. 





Date 


N 


T(s) 


IQ (") 


Airmass 


FORS1 


2000-11-21 


3 


3600 


0.74 (0.08) 


1.14 


FORS2 


2006-11-25 


5 


2915 


0.71 (0.10) 


1.24 




2007-01-16 


10 


5830 


0.93 (0.18) 


1.18 




2007-01-22 


5 


2915 


0.82 (0.11) 


1.19 




2007-02-11 


5 


2915 


0.82 (0.12) 


1.30 



January 16th and 22nd, and February 11th 2007 (see Tab. 
1 for a summary). The 2000 observations were performed 
with FOcal Reducer/low dispersion Spectrograph (FORS1), 
a multi-mode camera for imaging and long-slit/multi-object 
spectroscopy, as part of the ESO guaranteed time pro- 
gramme. At the epoch of the observations FORS1 was 
equipped with the original four port 2048x2084 CCD de- 
tector and it was mounted at the VLT Antu telescope. The 
observations were performed in standard resolution mode, 
with a 0'.'2 pixel size and a field of view of 6f 8 x 6f8. The 
low gain, fast read-out, single port mode was chosen. A se- 
quence of three 1200 s exposures was obtained through the 
Bessel B filter, with an airmass of ~ 1.14, an image qual- 
ity of ~ C/7, and dark time conditions. Since the seeing 
values measured by the differential image motion monitor 
(DIMM) arc relative to the zenith and not to the pointing 
direction of the telescope they are not necessarily indica- 
tive of the actual image quality. We thus computed the 
actual image quality from the measured point spread func- 
tion (PSF), derived by fitting the full width half maximum 
(FWHM) of a number of well-suited field stars using the 
S extractor tool (Bertin & Arnouts 1996), as documented in 
the FORS1 data quality control pagefl 
Two 600 s additional exposures were obtained in the V fil- 
ter but since the proposed candidate counterpart was not 
detected we focus our analysis on the B-band data only. 
Sky conditions were reported to be photometric (see Haberl 
et al. 2004 for a more detailed observations description). A 
very bright star, CD-50 1353 (B = 9.9, as listed in Simbad), 
located at « 45" from the position of RX J0420.0-5022 
was partially masked using the FORS1 occulting bars. Bias, 
twilight flat-fields frames, and images of the standard star 
fields SA 92 and Rubin 149 (Landolt 1992) were obtained 
as part of the FORS1 science calibration plan. 

The 2006/2007 observations were performed with 
FORS2 as part of the ESO open time programme. At 
the epoch of the observations, FORS2 had swapped with 
FORS1 at the VLT Antu telescope. FORS2 is equipped 
with two 2kx4k MIT CCD detectors. Due to vignetting, the 
effective sky coverage of the two detectors is smaller than 
the projected detector field of view, and it is larger for the 
upper CCD chip. Observations were performed in high res- 
olution mode, with a 2x2 binning and a pixel size of 0'.'125. 
The low gain, fast read-out mode was chosen. The telescope 
pointing was set in order to position RX J0420. 0—5022 in 

2 http:/ /www. eso.org/observing/dfo/quality/FORSl/qc/qcl. html 
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the upper CCD chip to include a larger number of reference 
stars for a precise image astrometry thanks to its larger ef- 
fective sky coverage (3.'5 x 2'). Sequences of 580 s exposures 
were obtained through the Bessel B filter. The bright star 
CD-50 1353 was more efficiently masked both by position- 
ing it at the centre of the gap between the two chips and 
by using the FORS2 occulting bars. Unfortunately, the dis- 
tance of RX J0420.0-5022 from the gap (~ 17") and the 
width of the occulting bars (~ 25"), which can move along 
one CCD direction only, made it impossible to mask com- 
pletely the star halo. Exposures were taken in dark time and 
under mostly clear but not perfectly photometric sky condi- 
tions. In particular, the night of February 11th was affected 
by the presence of thin variable cirri. Atmospheric condi- 
tions were not optimal either. The first two nights were 
affected by a strong wind, close to the telescope pointing 
limit, while the nights of January 16th and February 11th 
were affected by 40 % humidity. Unfortunately, although 
foreseen by the instrument science calibration plan, only 
for some nights both day and night time calibration frames 
were taken. In particular, no twilight flat -fields were taken 
for the night of November 25th, while B-band standard star 
images (of the Rubin 152 field) were taken on the night of 
February 11th only. 

2.2. Data reduction and calibration 

We retrieved the FORS1 science images from the pub- 
lic ESO archival and we reduced them using tools avail- 
able in MIDAS for bias subtraction, and flat-field correc- 
tion. The same reduction steps were applied to the FORS2 
science images through the ESO FORS2 data reduction 
pipelineQ. We searched the archive for suitable twilight flat- 
fields to reduce the FORS2 November 25th science im- 
ages but the closest in time were those associated with our 
January 16th images. We evaluated the possibility of us- 
ing lamp flat-fields as backup calibration frames, with the 
caveat that they arc affected by reflections produced by 
the instrument atmospheric dispersion correctors. However, 
since lamp flat-fields are only taken for trending purposes, 
and the high resolution mode is not the standard one for 
FORS2, no suitable data was found in the archive. Thus, 
since the November 25th science images can not be cali- 
brated with twilight flat-fields taken on the same night, 
initially we do not use them in the subsequent analysis. 
For both the FORS1 and FORS2 data sets the photomet- 
ric calibration was applied using the available, extinction 
corrected, night zero points available through the instru- 
ment data quality control databas^E For the January 16th 
and 22nd FORS2 observations, for which no standard star 
images were taken, we assumed as a zero point the value 
extrapolated from the night zero point trend. Since none 
of the FORS2 observations was taken in perfectly photo- 
metric conditions, we estimated that a relative photome- 
try calibration, with the February 11th observations taken 
as a reference, would introduce an uncertainty comparable 
to that associated to the extrapolation of the zero point 
trend. We converted the trended FORS1 and FORS2 zero 
points, computed in units of electrons/s, to units of ADU/s 
by applying the corresponding electrons-to-ADU conver- 

3 http: / / archive.eso.org 

4 http:/ /www. eso.org/observing/dfo/quality/FORS2/pipeline 

5 http:/ /www.eso.org/observing/dfo/quality/FORS2/qc/qcl.l 



sion factors. For each of the two data sets, we then used 
the MIDAS task average/window to cosmic rays filter and 
stack single exposures. 

2.3. Astrometry 

As a reference for the astrometric calibration we used the 
GSC-2 version 2.3 (Lasker et al. 2008). Approximately 70 
GSC-2 objects are identified in the 2000 FORS1 image. 
From this list we filtered out extended objects, stars that 
are either saturated or too faint to be used as reliable astro- 
metric calibrators or too close to the CCD edges. We finally 
performed our astrometric calibration using 30 well-suited 
GSC-2 reference stars, evenly distributed in the FORS1 
field of view. The pixel coordinates of the selected GSC- 
2 stars were measured by fitting their intensity profiles 
with a Gaussian function using the dedicated tool of the 
Graphical Astronomy and Image Analysis (GAIA) inter- 
faceQ. The coordinate transformation between the detector 
and the celestial reference frame was then computed using 
the Starlink package ASTROI"0 using higher order polynomi- 
als to accounts for the CCD distortions. The rms of the 
astrometric solution turned out to be « 0?.'2, accounting 
for the rms of the fit in the right ascension and declina- 
tion components. Following Lattanzi et al. (1997), we esti- 
mated the overall uncertainty of our astrometry by adding 
in quadrature the rms of the astrometric fit and the preci- 
sion with which we can register our field on the GSC-2 ref- 
erence frame. This is estimated as y/3x o-gsc/VN s > where 
the y3 term accounts for the free parameters (x-scale, y- 
scale, and rotation angle) in the astrometric fit, ugsc is 
the mean positional error of the GSC-2 coordinates (0'.'3, 
Lasker et al. 2008) and N 8 is the number of stars used for 
the astrometric calibration. The uncertainty on the refer- 
ence stars centroids is below O'.'Ol and was neglected. We 
also added in quadrature the O'.'IS uncertainty (Lasker et al. 
2008) on the tie of the GSC-2 to the International Celestial 
Reference Frame (ICRF). Thus, the overall accuracy of the 
FORS1 astrometry is 0'.'27 (ler). The astrometric calibra- 
tion of the FORS2 image was computed in the same way but 
with a lower number of reference stars due to the smaller 
field of view of the FORS2 chip. The rms of the astrometric 
fit then turned out to be 0'.'36. Again, after accounting for 
systematic uncertainties (see above) the overall accuracy of 
the FORS2 astrometry is 0"43 (la). 

As a reference to compute the RX J0420. 0—5022 posi- 
tion we considered X-ray coordinates derived from Chandra 
observations which are closest in time to our VLT ob- 
servations. In particular, for the FORS1 observations 
(epoch 2000.89) we used the Chandra coordinates (epoch 
2002.86) published in Haberl et al. (2004), i.e. aj 20 oo = 
04 fe 20 m 01.95 s , 5j2ooo = -50°22'48'/l which have a nom- 
inal error of 0"6 (90% confidence level). For the FORS2 
observations (epoch 2007.04) we reanalysed a more recent 
Chandra observations (epoch 2005.85). Like in Haberl et al. 
(2004), we determined the source position with the CI AO 
task celldetect and we obtained aj20oo = 04 ?l 20 m 01.94 s , 
<5j2ooo = -50°22'48'.'2 (0"6; 90% confidence level). As 
shown in Haberl et al. (2004) a match between the coordi- 
nates of the X-ray sources detected in the Chandra field 
with those of their possible USNOB.10 counterpart did 

6 star-www.dur.ac.uk/ pdraper/gaia/gaia.html 
nl 7 http:/ /star- www. rl.ac.uk/Software/software. htm 



4 



Mignani et al.: Optical observations of RX J0420.0-5022 




Fig. 1. Left panel: 1' x 1' cutout of the coadded F0RS1 B-band image (3600 s) of the RX J0420.0-5022 field. North 
to the top, east to the left. The RX J0420. 0—5022 position lies at the centre of the 10" x 10" square and is indicated 
by the circle. The radius of the circle (0"87; 90% confidence level) represents the uncertainty on the computed 2000 
RX J0420. 0—5022 position (see Sect. 2.3), which accounts for the absolute accuracy of the Chandra coordinates at the 
reference epoch, the uncertainty on the proper motion extrapolation at the observing epoch, and the accuracy of our 
astrometric calibration. The white triangle at the bottom of the image is the edge of the FORS1 occulting bar. Middle 
panel: Enlargement of the 10" x 10" region after sky background subtraction. The intensity scale has been adjusted for 
a better visualisation of the faintest objects in the field. Objects labelling is as in Haberl et al. (2004). The faint feature 
labelled / is their proposed candidate counterpart. Right panel: image of the same region smoothed with a Gaussian 
filter over cells of 3 x 3 pixels. 



not reveal any significant systematic shift. No significant 
shift is found between the coordinates of the same X-ray 
sources between the 2002 and 2007 Chandra observations 
either. Thus, no boresight correction was applied to our 
reference coordinates. Recently, an upper limit on the RX 
J0420. 0—5022 proper motion (123 mas yr^ 1 , 2a) was ob- 
tained with Chandra (Motch et al. 2009) . We accounted for 
the proper motion uncertainty when we registered the refer- 
ence Chandra coordinates on the FORS1 and on the FORS2 
images. This yields an additional position uncertainty due 
to the unknown proper motion of ~ 0"123 and ~ 0"073 
(1<t) for each of the two images, respectively. The overall 
uncertainty to be attached to the RX J0420. 0—5022 posi- 
tion at the epoch of the FORS1 and FORS2 observations 
was finally obtained by adding in quadrature the error on 
the Chandra coordinates (Iff), the coordinate uncertainty 
due to the proper motion, and the overall error of the as- 
trometric calibration. This yields to uncertainties of 0"87 
and (90% confidence level) on the RX J0420.0-5022 
position on the FORS1 and on the FORS2 image, respec- 
tively. 

3. Results 

3.1. The FORS1 observations 

We first re-analyzed the VLT observations taken in 2000 
to better assess the confidence of the optical identifica- 
tion of RX J0420.0-5022 proposed in Haberl et al. (2004). 
Fig. 1 (left) shows the computed Chandra position of RX 
J0420.0-5022 overlaied on a cutout of the 2000 FORS1 
co-added B-band image. As seen from Fig. 1 (left), only 
the bright PSF core of star CD-50 1353 is masked, while 
its halo extends close to the target position. This increases 
the local sky background as well as the background noise, 
which results in a larger number of spurious detections. In 
order to enhance the detection significance for fainter ob- 



jects we tried to minimise the effects of the halo of star 
CD-50 1353 on the local sky background. Firstly, we fitted 
the sky background in an area of ~ 15" x 15" around the 
target position using a second order polynomial and we sub- 
tracted the fitted value from the co-added B-band image 
using the MIDAS task f it/f lat_sky. We warn here that 
the fit to the sky background is biased by the choice of the 
sampling areas. This can yield to more or less evident fea- 
ture enhancements when the sky background subtraction is 
applied to the image. We thus carefully choose the sampling 
areas not to introduce systematic effects in our procedure. 
Fig. 1 (middle) shows a zoom of the sky-subtracted im- 
age. As already shown by Haberl et al. (2004), four objects 
are clearly detected close to the Chandra position. In addi- 
tion, a very faint feature is possibly recognised within the 
Chandra error circle. We identify this feature with object 
/ of Haberl et al. (2004), which they tentatively proposed 
as a candidate counterpart to RX J0420. 0—5022. However, 
the excess of counts at the feature position is comparable to 
the rms of the local sky background, which corresponds to 
a very low detection significance of rs 2ff. We re-computed 
the magnitude of the feature through PSF photometry. We 
derived the PSF parameters from a set of several non satu- 
rated objects selected for their stellar like profiles, located 
close to the RX J0420. 0—5022, and spanning a large range 
of magnitudes. The airmass correction was applied using 
the Paranal extinction coefficients measured with FORS^j. 
We found B = 27.52 ± 0.61. This is fainter than the value 
of B = 26.57 ± 0.30 reported in Haberl et al. (2004) but 
it is still compatible at the Iff level when systematic un- 
certainties in their photometry are taken into account (see 
section 3 of Haberl et al.). For a better visualisation, we 
smoothed the image using a Gaussian filter over cells of 
3x3 pixels i.e. of size comparable to that of the image PSF. 
Since the image smoothing enhances the detection of very 

8 http:/ /www. eso.org/observing/dfo/quality/FORSl/qc/qcl. html 
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Fig. 2. Left panel: 1' x 1' cutout of the co-added FORS2 B-band image (7540 s) obtained by the co-addition of the 
January/February 2007 best image quality exposures. North to the top, east to the left. The circle (l'.'ll radius; 90 % 
confidence level) corresponds to the uncertainty on the computed RX J0420. 0—5022 position (see Sect. 2.3). The white 
band in the left panel corresponds to the gap between the two CCD chips. Middle panel: zoomed image of the inner 
10" x 10" region (marked by a square) after co-addition of all the available exposures (see text), rebinning, and sky 
background subtraction. Right panel: image of the same region smoothed with a Gaussian filter over cells of 5 x 5 pixels. 
The faint feature detected within the Chandra error circle is labelled X. 



faint objects but also that of fluctuations of the noisy sky 
background we fine-tuned the smoothing parameters not to 
produce an over-enhancement of background features. The 
result is shown in Fig. 1 (right). However, the image pro- 
cessing (sky subtraction and smoothing) does not single out 
object / against the many, similarly significant, background 
features recognised around the Chandra position. Thus, we 
can not rule out that object / was a spurious detection due 
to the high background noise induced by the halo of star 
CD-50 1353. 



3.2. The FORS2 observations 

We used our follow-up FORS2 observations to search 
for a higher confidence candidate counterpart to RX 
J0420. 0—5022. In order to minimise the effects of the halo 
of star CD-50 1353, we first co-added only the exposures 
taken with an image quality better than 1". In first 
place, we used the co-addition of the best image quality 
exposures of all nights with the exception of those taken on 
November 25th 2006, which were calibrated using twilight 
flat -fields taken about 40 days apart (see Sect 2.2). 
Fig. 2 (left) shows the computed Chandra position of RX 
J0420.0-5022 overlaid on a cutout of the FORS2 B-band 
image (6960 s) obtained from the co-addition of the twelve 
best image quality (0"8-0"9) January/February 2007 expo- 
sures. Indeed, although the more efficient masking reduced 
the contamination from the halo of star CD-50 1353, 
the sky background at the Chandra position remained 
significantly affected by scattered light. As we did in Sect. 
3.1, we fitted and subtracted the sky background from 
the co-added image. In order to increase the S/N ratio 
per pixel we then rebinned the sky-subtracted image by a 
factor of 2, ending up with a pixel size of 0"25 which well 
matches that of FORS1 (0"2). We did not find evidence for 
object /, the feature tentatively proposed by Haberl et al. 
(2004) as a candidate counterpart to RX J0420. 0-5022. 
However, we possibly recognised a second feature within 
the Chandra error circle, « 0"5 north of the expected 



position of object /. Unfortunately, the low number of 
counts only yields to a marginal detection significance 
(~ 3a). As a test, and being aware of possible issues 
related to the non optimal flat -fielding, we decided to 
use the five November 2006 exposures (2915 s) which 
happen to have the best image quality ((f.'7) in the FORS2 
data set. As done for the January/February data set, 
we fitted and subtracted the sky background from the 
co-added image and we rebinned the sky-subtracted image 
by a factor of 2. Interestingly, a feature appears right at 
the same position of that seen in the co-addition of the 
January/February best image quality exposures, although 
with only a ~ 2.5tr detection significance. While we do not 
claim that this is a strong detection evidence, it is quite 
unusual that a background feature appears at the same 
position in images taken weeks apart. To increase the S/N 
ratio, we both co-added all the twenty 2007 exposures 
(11680 s) and all the available exposures (14575 s), again 
applying sky-subtraction and rebinning, and we obtained 
a detection significance of ~ 3.5<r and « 3.9cr, respectively. 
A zoom of the longest integration time, co-added image 
is shown in Fig. 2 (middle), where the feature detected 
in the Chandra error circle is labelled X. For a better 
visualisation, we smoothed the image using a Gaussian 
filter over cells of 5 x 5 pixels (Fig. 2, right). 

As done in Sect. 3.1, we measured the flux of object 
X through PSF photometry. The airmass correction was 
applied using the Paranal extinction coefficients measured 
with FORS^. Due to the still low S/N in the aperture 
and to the noisy sky background the flux measurement 
is obviously affected by a large error. Our best estimate 
gives B = 27.5 ± 0.3, where the statistical error obviously 
dominates over the uncertainty of our absolute photometry 
(Sect. 2.2). The object magnitude and coordinates are listed 
in Table 2 together with those of the other objects identi- 
fied in Fig. 2, as a reference. Interestingly, the flux of object 
X coincides with that of object I (B — 27.52 ± 0.61), the 

9 http:/ /www. eso.org/observing/dfo/quality/FORS2/qc/qcl. html 
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Table 2. Label, coordinates, and B-band magnitudes of 
the objects identified in the F0RS2 image (Fig. 2; middle). 
Coordinate uncertains are derived from our astrometric cal- 
ibration (Sect. 2.3). A photometry calibration error of 0.05 
magnitudes is assumed (Sect 2.2). 



ID 


(hms) 
"J2000 


U J2000 


B 


X 


04 20 01.94 


-50 22 47.75 


27.5 ± 0.3 


A 


04 20 02.10 


-50 22 42.60 


24.35±0.05 


D 


04 20 01.59 


-50 22 43.17 


24.95±0.05 


E 


04 20 01.57 


-50 22 48.75 


25.37±0.07 


H 


04 20 02.28 


-50 22 46.30 


26.49±0.10 



Blackbody 




10- 3 0.01 0.1 1 

Energy (keV) 



Fig. 3. Best fit models to the XMM-Newton EPIC-pn spec- 
trum of RX J0420. 0-5022 (green data points). The red and 
black lines correspond to a single blackbody and to a black- 
body plus an absorption line spectral fit, respectively (see 
Sect. 4 for details). Absorption-corrected model curves are 
drawn as dashed lines. The dereddened B-band flux of the 
putative candidate counterpart (la error) is marked. 

candidate counterpart tentatively proposed by Haberl et al. 
(2004), which we re-computed in Sect. 3.1. One may thus 
speculate whether we detected the same feature both in the 
FORS2 and in the FORS1 images, although at slightly dif- 
ferent positions. The measured angular separation between 
object X and object / is 0'.'5 ± 0'.'3, accounting for an es- 
timated uncertainty of one pixel on the object centroid in 
both the FORS1 and FORS2 images. This would imply a 
yearly displacement of 80 ± 50 mas yr _1 , consistent with 
the upper limit on the RX J0420. 0—5022 proper motion 
(Motch et al. 2009). The yearly displacement would thus 
imply a transverse velocity of sa 140 e?35o km s _1 , where 
d-350 is the neutron star distance in units of 350 pc (Posselt 
et al. 2007), i.e. within the range of the tangential veloc- 
ities inferred for neutron stars. The actual proper motion 
measurement, to be eventually obtained with Chandra, will 
unambiguously address this speculation. 

4. Discussion 

The very marginal detection significance (~ 3.9c) of ob- 
ject X against the number of local spurious detections 
makes it difficult to determine whether or not it is real and, 
thus, whether or not we have detected a candidate optical 
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Fig. 4. Upper panel: blackbody temperature T Q as a func- 
tion of the emission radius r Q for different values of the 
optical excess /. The red, blue, and green lines correspond 
to an optical excess of / = 7, 4 and 2, respectively. For each 
value of /, the three curves are drawn for different values of 
the neutron star distance, 550, 350 and 200 pc (dashed, dot- 
dashed, and solid lines, respectively). Lower panel: relative 
contribution R to the total 0.1-lkeV X-ray flux of a black- 
body with temperature T a versus the radius of the emitting 
region r Q for different values of the optical excess / and of 
the source distance. The line style and colour coding is the 
same as in the upper panel. The horizontal dotted line cor- 
responds to the threshold R = 0.1. The allowed region in 
the parameters space lies below this line. 



counterpart to RX J0420. 0-5022. We compared the flux 
of the putative candidate counterpart with the extrapola- 
tion in the optical domain of the models which best fit the 
XMM-Newton EPIC-pn spectrum of RX J0420.0-5022. To 
this aim, we have re-analysed the original data of Haberl 
et al. (2004) using updated calibration files. The spec- 
trum can be fit by a single blackbody with temperature 
kT x = 46.2 ± 1.4 eV and N H = (0.73 ± 0.21) x 10 20 cm" 2 
(reduced x 2 = 2.08, 64 d.o.f.), corresponding to an emission 
radius rx = 5.11 ^350 km, where rx is the X-ray emission 
as seen from infinity and ^350 is the neutron star distance 
in units of 350 pc (Posselt et al. 2007). However, a black- 
body with kT x = 47.8 ± 2.2 eV and N H =1.19tg;|i x 10 20 
cm -2 (rx = 5.44 g?35o km) plus an absorption line with 
centroid energy Eu ne = 337 ± 24 eV and equivalent width 
EWu n e = 47±5 eV gives a better fit (reduced x 2 = 1-33, 62 
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d.o.f.). The line width au ne ) was fixed at 70 eV, as in Haberl 
et al. (2004). For the spectral fits we used element abun- 
dances both from Anders & Grevesse (1989) and Wilms et 
al. (2000), obtaining virtually the same results. The best- 
fit, absorption-corrected X-ray spectra of RX J0420. 0—5022 
are shown in Fig. 3 together with the optical flux of its pu- 
tative counterpart. We corrected for the absorption in the 
B band using as a reference the Nh derived from the best-fit 
X-ray spectral model (blackbody plus absorption line) and 
applying the relation of Predehl & Schmitt (1995) with the 
extinction coefficients of Fitzpatrick (1999). From the flux 
of the putative counterpart we can rule out a > 7 optical ex- 
cess with respect to the extrapolation of the XMM-Newton 
spectrum. We note that an optical excess of ~ 5 is usually 
observed in other optically identified XDINSs with the ex- 
ception of RX J1605.3+3249 and 1RXS J214303.7+065419, 
where it is as large as ~ 15 (Motch et al. 2005) and w 30—40 
(Zane et al. 2008; Schwope et al. 2009), respectively. 
As a limit case, we checked whether an optical excess of ~ 7 
would be compatible with either rotation-powered emission 
from the neutron star magnetosphere or with thermal emis- 
sion from a fraction of the neutron star surface, colder and 
larger than that responsible for the X-ray emission. In the 
first case, the value of the X-ray period and the upper limit 
on the period derivative of RX J0420.0-5022 (P = 3.45s; 
P < 92 x 10~ 13 s s" 1 ; see Haberl 2007) only yield a ro- 
tational energy loss E < 8.8 x 10 33 erg s _1 . The flux of 
the putative counterpart would imply an optical luminosity 
Lb ~ 1.2 x 10 27 erg s _1 d% 50 . This would correspond to an 
emission efficiency t]b = Lb/E > 1.3 x 10~ 7 , which could 
still be compatible with the values expected for 10 6 — 10 7 
years old neutron stars (Zharikov et al. 2006). However, a 
period derivative P ~ 10~ 13 s s _1 , comparable to that of 
other XDINSs, would imply a factor of 100 lower E and 
would make it less likely that the optical emission is pow- 
ered by the rotational energy loss. In the second case, we 
can constrain both the blackbody temperature T Q and the 
emission radius r a , as seen from infinity. Since the fit to the 
XMM-Newton spectrum does not require the presence of 
a second blackbody component at lower temperature, we 
can impose that its relative contribution R to the total X- 
ray flux in the 0.1-1 keV band (see Sect. 3 of Zane et al. 
2008) must be << 1. We chose R = 0.1 as a reasonable 
threshold. We first computed the values of T Q for a grid 
of values of r Q and for different values of the optical ex- 
cess / = ^ y° an d °f the source distance (Fig. 4, upper 

panel), where Tx = 47.8 cV and rx = 5.44 ^350 km are 
derived from the best X-ray spectral fit (blackbody plus 
absoption line, see above). We then computed R from the 
values of r a and T Q (Fig. 4, lower panel). As it is seen, for 
a neutron star distance of 350 pc an optical excess of ~ 7 
would be compatible with a blackbody with kT a < 25 eV 
and an an implausibly large emitting radius of r Q > 23 km. 
Thus, would our putative counterpart be confirmed, an op- 
tical excess of ~ 7, for a neutron star distance of ~ 350 
pc, might rather point towards a non-thermal origin for 
the optical emission, as proposed for RBS1774 (Zane et al. 
2008). Actually, as Fig. 4 shows, unphysical large radii are 
required even if the actual counterpart is dimmer, / ~ 2-4, 
unless the neutron star is at ^ 200 pc. 



5. Conclusions 

We carefully re-analysed archival VLT/ FORS1 observa- 
tions of the field of the XDINS RX J0420.0-5022, taken 
in 2000, and we performed deeper follow-up observations 
with FORS2 in 2006 and in 2007. With a measured de- 
tection significance of ~ 2a and a re-computed flux of 
B = 27.52 ± 0.61, we can not rule out that the candi- 
date counterpart tentatively detected in the FORS1 im- 
ages by Haberl et al. (2004) was a feature of the noisy 
sky background, produced by the very bright nearby star 
CD-50 1353. While we could not confirm this detection in 
our deeper FORS2 images, we detected an apparently new 
feature (B = 27.5 ± 0.3) within the updated Chandra error 
circle of RX J0420.0-5022, w Of! 5 north of the expected po- 
sition of that detected in the FORS1 images. Interestingly, 
both their similar flux and their angular separation, com- 
patible with the upper limit on the RX J0420.0-5022 
proper motion, suggest that we might have actually de- 
tected the same feature both in the FORS1 and in the 
FORS2 images. However, its still marginal detection sig- 
nificance (~ 3.9ct) makes it difficult to determine whether 
the latter feature is associated with a real object, and thus 
it is the RX J0420. 0—5022 candidate counterpart, or it is 
also a possible background feature. From the flux of the pu- 
tative counterpart we can rule out a > 7 optical excess with 
respect to the extrapolation of the XMM-Newton spectrum. 
An optical excess of <~ 7 (or lower) could be compatible ei- 
ther with rotation-powered emission from the neutron star 
magnetosphere or with thermal emission from the neutron 
star surface for a distance 200 pc , i.e. much lower than 
the current best estimate of <~ 350 pc (Posselt et al. 2007). 
More observations are required to confirm the detection of 
the putative candidate counterpart and to validate its iden- 
tification with RX J0420. 0-5022. Unfortunately, the pres- 
ence of star CD-50 1353 severely hampers ground-based 
follow-up observations, even if performed under sub-arcsec 
seeing conditions and using a very careful masking. High 
spatial resolution observations with the refurbished Hubble 
Space Telescope (HST), possibly to be performed in the ul- 
traviolet, are the only way to settle the identification issue. 
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